Inhomogeneous dynamic nuclear polarization and suppression of electron-polarization 

decay in a quantum dot 
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We investigate the dynamic nuclear polarization process by frequently injecting polarized electron 
spins into a quantum dot. Due to the suppression of the direct dipolar and indirect electron-mediated 
nuclear spin interactions, by the frequently injected electron spins, the analytical predictions under 
the independent spin approximation agree well with quantum numerical simulations. Our results 
show that the acquired nuclear polarization is highly inhomogeneous, proportional to the square of 
the local electron-nuclear hyperfine interaction constant, if the injection frequency is high. Utilizing 
the inhomogeneously polarized nuclear spins as an initial state, we further show that the electron- 
polarization decay time can be extended 100 times even at a relatively low nuclear polarization 
(< 20%), without much suppression of the fluctuation of the Overhauser field. Our results lay the 
foundation for future investigations of the effect of DNP in more complex spin systems, such as 
double quantum dots and nitrogen vacancy centers in diamonds. 

PACS numbers: 03.65.Yz, 03.67.Pp, 73. 21. La 



I. INTRODUCTION 

Quantum coherence is of key importance for practi- 
cal quantum devices, such as quantum computerai*^ and 
spintronic devices^. To maintain the coherence of a quan- 
tum system, there are basically two kinds of methods: 
(I) Effectively reducing the system and environment cou- 
pling strength such as dynamical decoupling^i^ and em- 
ploying decoherence free subspace^ii; (II) Engineering 
the environment to a certain states, e.g., a dark state, 
a many-body singlets, which only weakly affect the sys- 
tem dynamics*^". 

Electron spins in semiconductor quantum dots (QDs) 
are promising candidate of solid-state quantum bits 
(qubits), due to their long relaxation time in strong 
magnetic fields^ii^ii^. However, the QD electron spins, 
which are decohered by surrounding nuclear spins 
through hyperfine coupling, have a rather short coher- 
ence/relaxation time, typically in the order of 10 ns in 
a few mT magnetic field and at '^100 mK low temper- 
atureiii^— . Thus, extension of electron spin coherence 
time is highly desired^Sr^^^. 

Many methods have been proposed to prolong the co- 
herence time of the electron spins, including dynamical 
decoupling via coherent control pulse o^^'^'^" — and quan- 
tum environment engineering (nuclear spin state narrow- 
ing and nuclear spin polarization)2r— iSiiS^i^^. For the 
dynamical decoupling and nuclear spin state narrowing, 
theoretical results show that the coherence time can be 
extended more than 2 orders of magnitude^ '^^i'^" , but the 
experimental confirmation is still illusive, due to the chal- 
lenging requirement of the high repetition rate of the con- 
trol pulses or the measurements. 

Nuclear spin polarization has been extensively inves- 
tigated both theoretically and experimentally. In gen- 
eral, the nuclear spins can be polarized by applying 
a very strong but static magnetic field (usually above 



1 Tesla)^ii^i^i2^ or by repeatedly injecting polarized 
electron spins with optical or electrical pumping meth- 
ods^>^i^~— . We refer the former as static nuclear po- 
larization and the latter as dynamic nuclear polarization 
(DNP). For the static nuclear polarization, the electron 
spin coherence time is almost unchanged in a QD with 
uniformly polarized nuclear spins, unless the nuclear po- 
larization is above 90% which is unexpectedly high in 
experimentsi^ii^'^i'^. For the DNP, early experiments 
and theories show that the coherence time of the elec- 
tron spin in double QDs can be extended significantly by 
suppressing the fluctuation of the Overhauser field dif- 
ferenceiSiSii^ilLj but subsequent experiments and more 
recent theories support a different interpretation^^—. A 
complete physical understanding of the DNP process, 
rooted in a microscopic picture, is still lacking. 

In this paper, we focus on the DNP process in a QD 
where the hyperfine coupling between the central elec- 
tron spin and the surrounding nuclear spins is domi- 
nant. We develop analytical equations for the DNP dy- 
namics under the independent spin approximation (ISA). 
Numerical quantum simulations with small-scale nuclear 
spin bath are further performed to confirm the validity 
of the analytical predictions. We ascribe the success of 
ISA to the suppression of the direct dipolar and indi- 
rect electron-mediated nuclear spin interactions by fre- 
quently injected electron spins. Furthermore, we find 
that the nuclear polarization acquired through DNP pro- 
cesses is highly inhomogeneous: strongly coupled nuclear 
spins gain much more polarization. Consequently, the 
electron-polarization decay is suppressed significantly, if 
the final DNP state is employed as an initial nuclear spin 
state, due to the protection effect by these highly po- 
larized nuclear spins, instead of the suppression of the 
Overhauser field fluctuation. 

The paper is organized as follows. We describe the QD 
system in Sec. |lTl The dynamics of a DNP process for a 
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single nuclear spin and many nuclear spins are presented 
in Sec. IIIII In Sees. llVlandlVlwe discuss respectively the 
effect of external magnetic fields and the dipolar cou- 
pling between nuclear spins during a DNP process. In 
Sec. IVII the suppression of the electron polarization de- 
cay is shown with the initial nuclear spin state being an 
inhomogeneously polarized one. Finally, the conclusion 
is given in Sec. IVIIl 

II. SPIN SYSTEM IN A QD 

We consider a singly charged quantum dot where the 
electron spin couples with the surrounding nuclear spins 
via Fermi contact hyperfine coupling 

N N 
k=l k=l i^k 

(1) 

where wo is the electron Zeeman splitting in an exter- 
nal magnetic field along z direction, S (S* = 1/2) is 
the electron spin, and Ik [Ik — 1/2, fc = 1,2,-- - 
with N the total number of nuclear spins) is the fcth 
nuclear spiniiii^i^. For an electron in a QD, Af. = 
(87r/3)(/|^B(7„^„|0(xfe)p is the Fermi contact hyperfine 
coupling constant, which is proportional to |(/)(xfc)|^, the 
electron density at the fcth nucleus. Here g* and gn are 
the Lande factors of the electron and the nuclei, respec- 
tively. The dipolar coupling strength between nuclear 
spins j and k are characterized by Tjk- We set fi = 1 for 
simplicity. The Zeeman terms of nuclear spins have been 
removed by adopting a rotating reference frame with the 
same frequency as that of the nuclear Larmor precession. 

For a gate-defined GaAs QD''^, the number of nuclear 
spins is about 10^ ^ 10^. The Zeeman splitting wq for 
the electron spin ranges from to 200 ^eV, proportional 
to the external magnetic field. Due to the inhomogeneity 
of the hyperfine coupling constant Ak is in general 
non- uniform and is typically in the order of Ak ^0.1 
neV. The dipolar coupling strength (nearest neighbor) is 
about Tjk ^ 0.01 neV, much smaller than Ak- We neglect 
the small nuclear dipolar interaction, i.e., Tjk = 0, unless 
otherwise stated. 



III. DNP 

The DNP process is consisted of many cycles. A typical 
DNP cycle can in principle be divided into approximately 
three steps as shown in Fig. [1] 

(I) Inject a polarized electron spin into the QD. The 
density matrix of the total system becomes 

Po=Po®Po, (2) 

where Pq= | t)(t I denotes the polarized electron 
spin state and denotes the nuclear spin state. 
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FIG. 1: (Color online) Schemetic of DNP process in a QD. 
S denotes the source of polarized electron spins (blue arrow), 
QD the quantum dot, and D the drain of the electron. (I) 
Injection of a polarized electron spin into the QD. (II) Due 
to the nonuniform hyperfine interaction between the electron 
and the nuclei in the QD, the electron polarization is inhomo- 
geneously transferred to the nuclear spins. (Ill) Ejection of 
the depolarized electron. After a DNP cycle, the total nuclear 
polarization is increased. The polarization of the strongly in- 
teracting nuclear spins with the electron increases faster than 
that of the weakly interacting ones. 



(II) Let the coupled spin system mix and the electron 
polarization is transferred to the nuclear spins, due 
to the total polarization conservation of the hyper- 
fine coupling. The evolved density matrix after a 
cycle time t is 

p{r) = lt)(tl®Pn + lt)ai®Pn 

Mi){n®ph + \\){i\®pli- (3) 

(III) Eject the electron out of the QD*^. The nuclear 
spin state after a DNP cycle then becomes 

P"(r) - (t|p(r)|t) + aip(r)U) 

= Pn + Pli- (4) 

Repeat the steps (I) — (III) till the nuclear spins reach 
a required polarization. In experiments^i^l, the DNP 
cycle period t can be as small as 250 ns, corresponding 
to a frequency of 4 MHz. 
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A. DNP of a nuclear spin 

As a starting point, we consider a single electron spin 
coupled to a nuclear spin. The Hamiltonian without a 
magnetic field is simplified as 



n 



AS ■ I. 



(5) 



The injected electron is fully polarized with a polarization 
Ps = '2{Sz) — 1- The initial nuclear polarization is p = 
'2{Iz)- The initial state of the total system is given by 



Po = |t)(tl® 



1 



It is easy to analytically find that the evolution of the 
electron polarization i&i^ 



1 



p 



1-p 



cos{At). 



(6) 



Due to the conservation of the total polarization, the 
polarization gained by the nuclear spin after a DNP cycle 
with a cycle period r is 



Ap = l-ps 
l-p 



[1 - cos(Ar)]. 



(7) 



If Ap <^ 1 for a single DNP cycle, or equivalently, 
At <C 1, the nuclear polarization difference Ap in the 
above equation can be approximated by a differential dp. 
We then obtain the following differential equation for the 
nuclear spin polarization during the DNP process 



dp 
dn 



A 



2^2 



(8) 



where n denotes the number of DNP cycle. Note that 
we have also taken the approximation cos(Ar) « 1 — 
since At ^ 1. 
Assuming the nuclear spin is initially unpolarized and 
solving the equation (|S]), we obtain the nuclear polariza- 
tion after n cycles is 



P{t 



1 



(9) 



It is clearly shown in Fig. [5] that the nuclear polar- 
ization exponentially approaches its saturation value of 
1. At the initial stage of the DNP process, where 
A^T'^n/A <^ 1, the nuclear polarization can be well ap- 
proximated by p « which indicates that the 
nuclear polarization is proportional to the square of the 
coupling strength. At the nearly saturated stage, where 
A^t'^u/A ^ 1, we may approximate the nuclear polariza- 
tion as p « tanh(A^T^n/8), as previously usedi^i^. We 
remark that the analytical prediction deviates from the 
exact numerical calculation once the condition At 1 
is seriously violated, as shown in the case of At — 1.0 in 
Fig. [2 
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FIG. 2: DNP of a nuclear spin with At = 0.2 (circles), 0.5 
(squares), and 1.0 (crosses). The solid lines are obtained from 
Eq. Q. 



If Ap ^ 1 is not satisfied, we get the following re- 
cursion relation between the nuclear polarizations of the 
(n — l)th cycle and the nth cycle 

Pn=Pn-l + - 1^ [1 - COs(At)]. (10) 



It is straightforward to find that 

'1 + cos(At) 



Pv 



1 - 



(11) 



where po = has been adopted. This is an exact solution. 
In the large n limit, the nuclear polarization becomes 

= 1 _e-[l-cos(Ar)]«/2^ (12) 

It is easy to obtain Eq. (0) by further assuming At <^ 1. 

B. DNP of a few nuclear spins 

It is challenging to exactly describe the DNP dynam- 
ics for many nuclear spins. By assuming the same cou- 
pling constant for all the nuclear spins, the DNP process 
in double quantum dots has been numerically investi- 
gated^M^iiiS,. While for different hyperfine coupling con- 
stant, it is only possible to numerically investigate the 
DNP process for a few nuclear spins N < 15. 

We consider a few nuclear spins {N = 4, 12) coupled 
to the electron spin with Gaussian distributed random 
coupling constant A^. The DNP process is investigated 
by numerically solving the von Neumann equation for the 
total density matrix 



ih 



dpjt) 
dt 



(13) 
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FIG. 3: DNP of iV = 4 nuclear spins with A^t = 0.2 (top) 
and 1.0 (bottom), where Am is the largest coupling constant. 
The solid lines are obtained from Eq. ((Oj. The agreement 
between numerical and analytical results shows that the ISA 
is valid if At is small. 
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FIG. 4: The same as in Fig. [3] except that the number of 
nuclear spins is A*' = 12. 

during the spin transfer time t £ (nr, [n + 1]t). The 
injection of the polarized electron spin at t ^ nr is 

p{t^nT)^\t){t\^Pn{nT). (14) 

The ejection at the end of the DNP cycle t — {n + 1)t is 
a partial trace over the electron spin 

p„(t) = Ti-sipit)} 

= {t\p{t)\t) + {i\p{t)\i)- (15) 

Since the interactions among the nuclear spins are 
much smaller than the hyperfine interaction between the 
electron and the nuclear spins, we may neglect the nu- 
clear spin interactions, including the direct dipolar in- 
teraction and the indirect electron-mediated-nuclear-spin 
interaction, to the leading order. Such an approximation 



is hereafter referred as independent spin approximation 
(ISA). Under the ISA, each nuclear spin is polarized in- 
dependently by the electron spin so that the nuclear po- 
larization obeys the analytical prediction of a single spin 
[see Eq. dH)]: pk{t = nr) = I - exp{-AlT^n/4). Ob- 
viously, the ISA is valid if the the effect of the nuclear 
dipolar interaction and the indirect electron-mediated- 
nuclear-spin interaction is small. 

During the DNP process, we quantitatively investigate 
the polarization of the fcth nuclear spin pk for different 
DNP cycle time r. The evolution of nuclear spin po- 
larization during DNP process, as well as the analytical 
predictions from Eq. are presented in Fig. [3] and 0] 
for = 4 and 12, respectively. For both A^s, we choose 
two typical values of r: (i) AmT = 0.2 is in the small r 
region where the ISA is valid; (ii) AmT = 1.0 is in the 
large r region where the ISA is violated, especially after 
many DNP cycles. 

By comparing the numerical results with the ISA pre- 
dictions in both Figs. [3] and 21 we find good agreement 
at small AmT {Am is the largest coupling constant), in- 
dicating the validity of the ISA for a few nuclear spins. 
The good agreement between the ISA predictions and the 
numerical results implies that fast DNP cycles effectively 
decouple the indirect electron-mediated-nuclear-spin in- 
teraction. This decoupling effect becomes more promi- 
nent if the DNP cycle period is small, similar to the 
quantum Zeno effecl^^Ti^. Given typical experimental 
conditions, where Ak ~ 0.1 neV and t ^ 50 ns^i, we 
obtain A^t ~ 0.001 which is safely in the valid region of 
ISA. While for large AmT, the ISA predictions deviate 
from the numerical results. The more the nuclear spins 
are, the larger the deviation is, as shown in Fig. [3] and SI 

Figures [3] and S] also show that the nuclear polariza- 
tion is inhomogeneous: strongly coupled nuclear spins 
are polarized much faster than the weakly coupled ones. 
The polarization of these strongly coupled spins reaches 
nearly 100% while the polarization of those weakly cou- 
pled ones is still very small. These results agree quanti- 
tatively with previous predictions of the two-region QD 
model^. 

With current experimental techniques, it is a big chal- 
lenge to measure directly the polarization of each nuclear 
spin, in particular, spatially inhomogeneous nuclear po- 
larization, during the DNP process. But the Overhauser 
field B induced by the nuclear spin polarization and the 
fluctuation of the Overhauser field AB are reachable^i*^. 
Therefore, we quantitatively compute these properties for 
different DNP cycle time r, as shown in Fig.[5]for TV = 12 
nuclear spins. 

In general, the Overhauser field is a vector, B = 

EJ5„e„ and its fluctuation along an a- axes is 
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AS, 



'(B2)-B2, where 



Br, 



It is easy to check that if the initial nuclear spin state 
is unpolarized, then during the whole DNP process 



AS, 



ihy) 

AB„ 



0. Therefore, 



B„ 



0, and 



— (1/2)a/^^ A1 , which are constant of 
motion. While along the z-direction, Bz and AB^ are 
both time dependent, as shown in Fig.[SJ Under the ISA, 
it is straightforward to calculate that (see Appendix |^ 



Bz = AkPk 



-Pi) 



(16) 
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FIG. 5: The reduced Overhauser field (top) and the reduced 
fiuctuation (bottom) during a DNP process for A'' = 12 nu- 
clear spins with AmT = 0.2. The solid lines are calculated 
under the ISA with Eqs. ©, (HB]), and The couphng 

constant Ak is the same as that in Fig. ID 



where pk — I — exp^—Alr^n/A) is given by Eq. dH). 

We present in Fig. [5] actually the reduced quantities, 
which are normalized to their maximal value 



B'z 
ABi 



2B, 



2ABz 



(17) 



Once again, we find from Fig. [5] that the analytical ISA 
results are pretty close to the numerical results, indicat- 
ing that the ISA is a good approximation for the param- 
eter of AmT = 0.2. For smaller A^t as in experiments, 
one would expect much better agreement between the 
analytical and the numerical results. 

The rapid increase of B'^ at the small number of DNP 
cycles in Fig. [5] implies that the strongly coupled nu- 
clear spins are firstly polarized. Under the ISA, the 
increase rate is proportional to A| and to the number 
of DNP cycles n. Once these strongly coupled nuclear 
spins are nearly fully polarized, the increase of B'^ be- 
comes slow, as shown at the large number of DNP cycles, 
due to the smallness of the hyperfine constant Ak of the 
weakly coupled nuclear spins. Quite differently, the fluc- 
tuation shows slow decrease when the Overhauser field is 
small but rapid decrease when the strongly coupled nu- 
clear spins are approaching the fully polarized state. The 
reason lies in that the flu ctuation of the fleld is roughly 
proportional to ^1 — pf. for these strongly coupled nu- 
clear spins, which is approximately a constant (zero de- 
crease rate) at small number of DNP cycles [n < 50 in 
the bottom panel of Fig. \5\ since pk is small. Once the 
strongly coupled nuclear spins are nearly fully polarized, 
i.e. pk ~ 1 for some nuclear spins, the fluctuation de- 
creases signiflcantly. However, we notice that the fluctu- 
ation of the Overhauser fleld is only suppressed around 
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FIG. 6: The reduced Overhauser field (top) and the reduced 
fiuctuation (bottom) during DNP process for N = 10^ nuclear 
spins with AmT = 0.02. The results are predictions under ISA 
with Eqs. H16|) . and p7p . The coupling constant A^ is 
randomly distributed with a Gaussian distribution function. 



1/3 of its initial value even at B'^ « 0.8, which implies 
that the DNP is not an efficient way to suppress the fluc- 
tuation in a QD. 



C. DNP in a QD 



For a gate-defined GaAs QDi 



the number 

of nuclear spins is in the order of 10^ and the hyperfine 
coupling is spatially inhomogeneous. Directly simulating 
the DNP process without approximation in such a huge 
system is far beyond the capability of current computers. 
Instead, we adopt the ISA which has been proved to be 
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FIG. 7: (Color online) The effect of external magnetic fields 
during the DNP process of a single nuclear spin. The fields are 
tJo = (black circles), ojq = 5 (blue asterisks), loq = 10 (red 
triangles), and uq — 15 (green crosses), in units of A. The 
parameters are At — 0.2. The solid lines are calculated from 
the analytical expression Eq. (|19|) . The nuclear polarization 
increases slower in a larger magnetic field. 




DNP cycle 

FIG. 8: (Color online) The same as Fig. [7] except the nuclear 
spin number is = 4. The black dashed lines are the nuclear 
polarizations of the 4 spins for luq = 0, the blue circles are for 
Ldo = 5 Am, and the blue solid lines are calculated under the 
ISA from Eq. (fT9)) . 



cycle is small. 



a good approximation at small number of nuclear spins 
and at small AmT. 

The distribution of the coupling constants Ak is as- 
sumed in a Gaussian form 

= ,18, 

where = 1 is the maximal coupling constant at the 
center, Xk,yk are random number in the region [-1, 1], 
and Zfc is random number in the region [-0.1, 0.1]. The 
position of the fcth nuclear spin is given by (xk,yk, z^)- 
We set (7_L = 0.2 and az = 0.1 x (t_l = 0.02, which de- 
scribes a pancake shape QD. 

We plot the analytical results under ISA in Fig. [6] for 
N = l(f and A,„t = 0.02. Similar to the case of TV = 12, 
the Ovcrhauser field increases quickly at small B'^ and 
becomes slow once B'^ is above 0.5. The Overhauser field 
fluctuation decreases quickly at the middle region of the 
DNP process, around n ^ 2 x 10*. Furthermore, the 
suppressed fluctuation is about 1/4 at B'^ being roughly 
0.9, which shows that the DNP method is not good at 
suppressing the Overhauser field fiuctuation. 



IV. EFFECT OF EXTERNAL MAGNETIC 
FIELDS 

When an external magnetic field is applied during the 
DNP process, the evolution of the nuclear polarizations 
is changed drastically. For a single nuclear spin, it is easy 
to obtain the following analytical expression after n DNP 
cycles, if the change of the nuclear polarization in a DNP 



p„ = 1 - e-(^Vf2')[l-cos(nr)]n/2^ j^^g^ 

where £7^ = A"^ + ujq. For a given r, the nuclear po- 
larization pn decreases almost inverse-quadratically with 
oscillations as the Zeeman splitting of the external mag- 
netic field Wo increases. If the DNP cycle time r is small 
enough, r <IC we find field-independent nuclear po- 
larization 

p„ « 1 - e-^'^'"/4. (20) 

We present in Fig. [7] the analytical prediction from 
Eq. (jl9p and numerical results, which agree well with 
each other for the given parameters. Such an agree- 
ment between the analytical and numerical results in- 
dicates the validity of the adopted approximation, i.e., 
the change of the nuclear polarization in a DNP cycle 
is small. Compared to the zero field result, the nuclear 
polarization increases slower with an external magnetic 
field applied. Such a polarization suppression becomes 
severer in larger magnetic fields. 

We shown in Fig. |8] the numerical results for = 4 
nuclear spins and the analytical results from Eq. 
under the ISA. The good agreement between them indi- 
cates that the ISA is a good approximation. Compared 
with the zero magnetic field results, we observe the po- 
larization suppression in a magnetic field. Furthermore, 
the suppression effect becomes more significant for the 
weakly coupled nuclear spins than that for the strongly 
coupled ones, due to the larger reduced magnetic field 
LOo/Ak- In this sense, we may deduce that the applica- 
tion of a magnetic field during the DNP process actually 
enhance the inhomogeneity of the nuclear polarizations. 
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FIG. 9: (Color online) Effect of nuclear dipolar interaction 
on the nuclear polarization during DNP processes. The cycle 
times are A^t = 0.02 (left) and AmT = 0.2 (right). The 
dashed lines are the analytical results under ISA and the blue 
circles are numerical results. For those weakly coupled nu- 
clear spins, spin diffusion induced by the dipolar coupling 
overwhelms the DNP induced by the hyperfine coupling. 



V. EFFECT OF NUCLEAR DIPOLAR 
INTERACTION 

The dipolar coupling between nuclear spins usually 
cause spin diffusion, i.e., the polarization are transferred 
from high polarization spins to low polarization spins. 
Since the nuclear polarization due to electron-nuclear hy- 
perfine coupling during DNP process is highly nonuni- 
form, including the nuclear dipolar coupling introduces a 
competing process which reverses or at least hinders the 
nuclear polarization inhomogeneity. 

We compute the DNP process of = 4 nuclear spins, 
including both the hyperfine coupling and the nuclear 
dipolar coupling. The parameters are generated ran- 
domly: Ai = 0.7059, A2 = 0.3009, ^3 = 0.0089, A4 = 
0.4014, = 5.673 x lO^^, T13 = 9.729 x lO^^ = 
1.332 X 10-2, = 6.490 x lO^^ = 1.734 x lO^^, 
= 3.909x10-2. The typical dipolar couphng strength 
is about one order of magnitude smaller than the typical 
hyperfine coupling strength. 

We present the numerical results in Fig.lH which shows 
clearly that the spin diffusion indeed occurs during the 
DNP process. Contrary to the usual free spin diffusion 
process, where the high polarizations decrease and the 
low polarizations increase with total nuclear polariza- 
tion conserved, the spin diffusion during DNP does not 
cause the high polarization decrease while the low polar- 
izations indeed increase with total nuclear polarization 
increasing. In other words, the dipolar interaction in- 
duced spin diffusion is more important for those weakly 
hyperfine-coupled nuclear spins whose polarization are 
mainly transferred from highly-polarized nuclear spins 
instead of the electron spin. Moreover, similarity in the 
left and right panel of Fig. IH] indicates that reducing the 
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FIG. 10: (Color online) Electron polarization decay time 
versus reduced Overhauser field for A'^ — 256 nuclear spins 
with Gaussian widths a± = 0.2A''xa (blue line with crosses), 
OAN±a (green line with circles), and 0.6A'^_ia (red line with 
squares) . The electron coherence time can be extended signifi- 
cantly even at moderate inhomogeneous nuclear polarizations 
(about B'^ = 0.3). Inset: Typical electron polarization decay 
for a± = 0.47Vxa at different /3 = 0, 5, 10, 20, 40, 80, 160, from 
bottom to top. The dashed horizontal line denotes Uz — 0.5. 



DNP cycle period r does not change significantly the fi- 
nal polarization distribution (inhomogeneity) , except the 
spin diffusion process is slowed (the total evolution time 
t = riT in the left panel of Fig. IH] is 10 times as that in 
the right panel). 



VI. SUPPRESSION OF ELECTRON 
POLARIZATION DECAY VIA DNP EFFECT 

DNP was originally proposed to extend the coherence 
time of the electron spin in the QD by suppressing the 
fluctuation of the nuclear-spin Overhauser field. The ex- 
tension of the coherence time of electron spins in QDs was 
already demonstrated in double QD experiments and in 
some theoretical work c-'^°'2^i29i43-46,52 ^ However, as we 

show in Fig. [S] and HI the fluctuation of the Overhauser 
fleld is not suppressed signiflcantly after a DNP process, 
even at pretty high nuclear polarization. We also notice 
that in the experiment done by Bluhm et alM the sup- 
pression of the fluctuation of the Overhauser field is not 
inversely proportional to the extension of the coherence 
time, i.e., two thirds suppression of the fiuctuation causes 
approximately seven times extension of Tj* (see Fig. 3 in 
Ref.'**). Therefore, other mechanisms, such as the pro- 
tection effect of the nuclear spins strongly coupled to the 
electron and highly polarized during DNP—, or the nar- 
rowing of the nuclear spin statea^'^'^i^'^, might also 
play an important role. 

We investigate the above conjecture by numerically 
simulating the electron spin polarization decay with an 



8 




FIG. 11: (Color online) The nuclear bath size effect on the 
electron polarization decay time. The electron density relative 
widths are the same for all three cases, a±/{N±a) = 0.4. To 
reach the same decay time ri/2, the required polarization of 
a large bath is smaller than that of a small bath. 



initial inhomogeneous nuclear polarization acquired dur- 
ing a previous DNP process. The coupling constants Ak 
is assumed in a Gaussian distribution form 



Ak = An cxp - 



(xk - xq)^ + {yk - yof 
2al 



(21) 



where Aq is the maximal coupling constant at the shifted 
center {xo,yo) = (0.1,0.27)a, {xk,yk) = (n^ - N±,ny - 
N±)a is the position of the fcth nuclear spin in a two- 
dimensional square lattice with a being the lattice con- 
stant and Hx^y G [l,2N_i_] the index of the fcth nuclear 
spin. By increasing the width cj^, we can adjust the 
Ak's from a sharp distribution to a flat one. 

The initial state is a product state of a fully polarized 
electron spin state and a partially and inhomogeneously 
polarized nuclear spin state. The electron spin is ini- 
tially polarized along the z direction. The nuclear spins 
are prepared through the DNP process as described in 
Sec. mil We adopt the previously obtained ISA expres- 
sion pk = 1 — exp(— 2/3A^) with /3 an adjustable param- 
eter. For f3Al <^1, Pk ~ '^I^Al, which manifests the fact 
that Pk is proportional to the square of the coupling con- 
stant. For pAf. ^ 1, saturates at 1, which agrees with 
the limiting cases of an infinite number of DNP cycles. 

The coupled electron-nuclear spin system evolves un- 
der the Hamiltonian described by Eq. ([1}, but excluding 
the nuclear dipolar coupling and setting zero the Zeeman 
splitting of the external magnetic field. The numerical in- 
tegration is performed with the P-representation density 
matrix method^'^'^. This method has been numeri- 
cally demonstrated as an efficient method for many-spin 
dynamics^i^ and shows good agreement with the exact 
method based on Chebyshev polynomial expansion of the 
evolution operator, in describing the electron polarization 
decay^. 



The electron polarization (25*2) is monitored during 
the evolution. With a polarized nuclear spin bath which 
generates a nonzero Overhauser field in average, the elec- 
tron polarization does not decay all the way down to zero. 
We thus subtract the nonvanishing value at very long 
time and normalize the decayed electron polarization as 
follows^ 



(5.)(t)-(5,)(oo) 
(5.)(0)-(5.)(oo) 



(22) 



where {Sz){t) = Tr{Szp(t)} with p(t) the density matrix 
of the coupled system at time t. 

Typical electron polarization decay processes are 
shown in the inset of Fig. (TU] We see an initial quick 
dip at short times and slow decay at long times. To 
characterize the decay of Uz , we define a half-decay time 
constant T1/2 at which Uz = 1/2. The relation between 
Ti /2 and the reduced Overhauser field B'^ is presented in 
the main panel of Fig. [TU]for various distribution widths 
of Ak- At small B'^, the half-decay time constant is al- 
most independent of B'^; while in a moderate region of 
B'^ ~ 0.3, the half-decay time increases rapidly and the 
growth rate becomes lager for a narrower distribution of 
Ak- This shows a big contrast to the relation between 
the field fluctuation AB'^ and B'^ in Fig. [5] and El where 
AB'^ changes slowly while B'^ increases. 

To investigate the bath size effect on the electron polar- 
ization decay, we present in Fig. [TT] the relation between 
Ti/2 and B'^ for N = 64, 256, and 1024. The resuhs 
clearly show that the half-decay time T1/2 increases more 
rapidly at moderate B'^ for a larger bath size. More im- 
portantly, the transition position of B'^ , where the growth 
rate of T1/2 changes abruptly, becomes smaller with the 
bath size increasing. Given this trend, the required B'^ 
is far below 0.1 for a real QD with N ^ 10^ nuclear 
spins if the electron polarization decay time is extended 
significantly. 

To understand the mechanism of the electron polariza- 
tion decay, we count the number of nuclear spins whose 
polarization is above a certain value pc in the case of 
N — 1024 and B'^ = 0.52 (the total nuclear polariza- 
tion is about 22%). By setting pc = 95%, the number is 
102; by setting pc = 90%, the number is 121. Note also 
2^1/2 — 0.7 at B'^ — with the unpolarized initial nuclear 
state and T1/2 — 70 at B'^ = 0.52. The electron polar- 
ization decay time is extended about 100 times, which 
roughly equals to the number of highly polarized (and 
strongly coupled with the electron) nuclear spins, indi- 
cating that these special spins protect the electron spin 
from decay^S. 

Another mechanism to explain the extension of the 
electron polarization decay time is the nuclear state nar- 
rowing^. We plot in Fig. [T^ the amplitude distribution 
of the density matrix elements before (top) and after 
(bottom) a DNP process for = 4 nuclear spins in 
its coupling basis, which is the eigen basis of the cou- 
pling Hamiltonian H — S ■ Aklk- Given a diagonal 
density matrix, the coupled electron-nuclear spin system 
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FIG. 12: (Color online) Amplitudes of the Density matrix 
elements of the coupled electron-nuclear spin system in cou- 
pling basis before (top) and after (bottom) the DNP process. 
The number of nuclear spins is N = 4. The density matrix is 
concentrated along the diagonal line, indicating a narrowing 
of the state distribution, through the DNP process. 



spin interaction. The application of an external mag- 
netic field does not change the overall DNP behaviors, 
except slowing the DNP process in proportion for all nu- 
clear spins. The inclusion of the nuclear dipolar interac- 
tion causes weakly coupled nuclear spins slowly polarized 
through nuclear spin diffusion. 

The inhomogeneously polarized nuclear spins after a 
DNP process are further utilized to extend the decay time 
of the electron polarization. Due to the protection effect 
of the highly polarized part of the nuclear spins, the elec- 
tron polarization decay time can be extended 100 times 
at a rather low nuclear polarization (about 20% for 1,000 
nuclear spins and smaller for more nuclear spins). 
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Appendix A: Derivation of Eq. (|16|) 

From the definition of Bz , we obtain 
Bz = (^Afe/fc.) 

k 

k 



would never change so that the decay time of the elec- 
tron polarization is infinity. In this sense, the degree of 
the concentration of the density matrix on the diagonal 
line characterizes the effect of the narrowing state. By 
comparing the top and the bottom panel of Fig. [121 '^^ 
clearly see a narrowing of the state distribution after a 
DNP process, indicating qualitatively the extension of 
the decay time. But more detailed work need to be done 
in the future to interpret quantitatively the extension of 
the electron polarization decay. 



VII. CONCLUSION 

In summary, we obtain an analytical expression for the 
inhomogeneous nuclear polarization acquired during a 
DNP process under the independent nuclear spin approx- 
imation in a singly charged QD, where the nuclear spins 
coupled with an electron spin through contact hyperfine 
interaction. Numerical simulations agree well with the 
analytical predictions in the limit of short DNP cycles, 
due to the suppression of the electron-mediated nuclear 



For the field fluctuation, we first calculate that 



k ] 



k 3 



\ E^'+EE^^^^w 
\ k k j^k J 



Then the field fluctuation becomes 
ABz = J{Bl)-Bl 



= i/E^^i^- 
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